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Abstract 

We expand the existing polaron response theory, expressed within the Mori-Zwanzig projection 
operator formahsm apphcable for the transfer of arbitrary energy and zero momentum, to the the 
case of finite momentum exchange. A general formula is then derived which can be used to calculate 
the response of a system to a probe that transfers both momentum and energy to the system. The 
main extension of the existing polaron response theory is the finite momentum exchange which 
was not needed until now since it is negligible for optical absorption. However, this formalism 
is needed to calculate the response of the polaronic system consisting of an impurity in a Bose- 
Einstein condensate to Bragg spectroscopy. We show that the well-known features that appear in 
the optical absorption of the solid state Frohlich polaron are also present in the Bragg response of 
the BEC-impurity polaron. The f-sum rule is written in a form suitable to provide an independent 
consistency test for our results. The effect of lifetime broadening on the BEC-impurity spectrum is 
examined. The results derived here are discussed in the framework of an experimental realization 
consisting of a lithium impurity in a sodium condensate. 
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I. INTRODUCTION 



The experimental realization of impurities in a Bose-Einstein condensate (BEC) [lH3| 
and the possibility to produce quantum degenerate atomic mixtures has resulted in an 
interest in the physics of impurities in a quantum gas. This has led to an investigation of 
the change of the properties of the bare impurities as a result of the interactions with the 

nn nn 

condensate (such as the effective interaction M, ^1 and the effective mass 9- 111) and a study 

nn 

of the self-trapping of the impurity |12l . Il3l | . Also the system of a spin down fermion in a 



sea of spin up fermions, usually indicated as the Fermi polaron, has been investigated which 
resulted in a good agreement between theory 

it was shown that when the Bogoliubov approximation is valid the system of an impurity in 
a BEC can be mapped to the Frohlich polaron system 18l-l20l| which shows that this system 



can be added to the list of condensed matter systems that can be imitated in the context of 
quantum gases {21I. The experimental realization of quantum degenerate atomic mixtures 



in an optical lattice 22|-|24j| has led to the theoretical study of the polaronic properties of 



impurities in a BEC where the influence of the optical lattice is felt only by the impurities 
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26| or by both the impurities and the condensate 27 1. 



In the context of solid state physics the Frohlich polaron consists of a charge carier 
electron, hole) interacting with the phonons in an ionic crystal or a polar semiconductor 
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29|. This Frohlich polaron Hamiltonian can not be diagonalised exactly and several ap- 



proximation methods have been developed for it. The most general way to study the system 
is through a variational principle within the path integral formalism which was developed 



by Feynman 



30|. The internal excitation spectrum of the Frohlich polaron was revealed 



through the calculation of the optical absorption by Devreese et al. [31j, based on the 
Feynman path integral method and a response formalism introduced in Ref. 32|. Recently 
the excitation spectrum of the Frohlich polaron Hamiltonian was also studied numerically 
with Diagrammatic Quantum Monte-Carlo numerical techniques and these results showed 
deviations in the large coupling optical absorption linewidth with the theoretical results o " 
3l| . see Ref. 33|. At weak and intermediate coupling the analytical theory of Ref. 31 1 
agrees well with the numerical simulations of Ref. [331] and this optical absorption spectrum 



was also studied experimentally 



34| . So far the strong coupling regime could not be probed 



experimentally since the largest attainable Frohlich polaron coupling constant in any solid 
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is not large enough. The large coupling behavior remains an important question since a 
better understanding of the intermediate and strong coupling regimes might be useful to 
elucidate the role of polarons and bipolarons in unconventional pairing mechanisms e.g. for 



high-temperature superconductivity 



35|. Furthermore it was shown that at intermediate 



coupling the introduction of an ad hoc lifetime r for the eigenstates of the polaron model 
system is an improvement of the theory which is called the extended memory function for- 
malism 



If the Bogoliubov approximation is applicable the Hamiltonian of an impurity in a BEC 
can be written as the sum of the mean field energy of the hornogeneous condensate and the 
Frohlich polaron Hamiltonian that describes the fluctuations 20|: 
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(1) 



which describes the interaction between the impurity and the Bogoliubov excitations where 
r and p represent the position and momentum of the impurity with mass m/, bj: and bt are 
the annihilation and creation operators for the Bogoliubov excitations, is the Bogoliubov 
frequency: 

ck\/l + m'/2, (2) 



and Vr is the interaction amplitude: 



1/4 



9lB- 



(3) 



In the above expressions use was made of the definition of the healing length of the conden- 
sate: C, = 1/ y/STcasBno with riQ = Nq/V the condensate density and asB the boson-boson 
scattering length. We also used the expression for the speed of sound in the condensate: 
c = h/ (^-s/2mB^)- In j2o| the Feynman variational path integral technique was applied to 
this Hamiltonian and an upper bound for the free energy was calculated by using the Jensen- 
Feynman inequality with the Feynman model system. This model system consists of the 
impurity mass coupled to another mass M by a spring with frequency fl. The parameters 
M and Q are then used to minimize the free energy. It followed that the Frohlich polaron 
coupling parameter in this case is given by: 



with aiB the impurity-boson scattering length. Together with the ratio between the masses of 
the impurity and the bosons this couphng parameter fully determines the static properties 
of the specific BEC-impurity polaron system. Depending on the value of this coupling 
strength two regimes were identified: a strong coupling regime which has properties that 
suggested a polaronic self-trapped state and a weak coupling regime which suggested a free 
polaron. Since this coupling parameter depends strongly on the scattering lengths, which 
can be tuned externally by a magnetic field through a Feshbach resonance (see eg. [37|), it 
may be possible to tune the system experimentally to the regime of strong coupling. This 
technique might enable the experimental realization of the strong coupling regime and reveal 
the internal structure of the Frohlich polaron Hamiltonian. 

Since, as opposed to the original Frohlich polaron in solid state, the BEC-impurity polaron 
system is not charged it is not helpful to perform optical absorption measurements on this 
system. The method to probe the system is addressed in the present paper and will be 
shown to be Bragg spectroscopy. This has proven to be a very succesfull tool to probe the 
structure of BEC's (see eg. 38| and jsQl). It is realized by probing the system with two laser 
beams with wave vectors ki and k2 and frequencies uji and uj2- An atom can then undergo a 
stimulated scattering event by absorbing a photon from beam 1 and emitting it in beam 2 
which changes its momentum hj k = ki — k2 and its energy hj = fiwi — huj2- See figure 
[T]for a schematical picture of a typical experimental set-up. A possible way to measure the 
response of the system is by performing a time-of-fiight experiment and count the number 
of atoms N Bragg that have gained^ a momentum k. Within the formalism of linear response 
theory this number is given by [37|]: 
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(5) 



with r the duration of the pulse, V the amplitude of the laser induced potential and x (^^, 
the density response function, defined by: 
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\pk\r^)\ 
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UJ — UOr. 



(6) 



with the induced density of the impurity, Z the partition function, \n > the eigenstates 
of the unperturbed system and ujnm = {E^ — Em) the transition frequencies. This density 
response function can be rewritten as the Fourier transform of the retarded density-density 



condensate 




FIG. 1: Schematical picture of a typical experimental set-up for Bragg spectroscopy. Two lasers 
with momentum ki and k2 and energy /kJi and fiuj2 are impinged upon the impurity in the con- 
densate which gains an energy fkjj = hui — hu}2 and a momentum k = ki — k2- 



correlation function: 



dte 



iujt 



(7) 



The retarded density-density correlation function is needed to obtain the response of a 
system to Bragg spectroscopy. 

We will start by using the Mori formalism to calculate a general expression for the density- 
density correlation function and then apply it to the generic polaron system with ar- 
bitrary dispersion and interaction amplitude. Then we show that if we calculate from this 
expression the current- current correlation for the k = limit we find the known result for the 



optical absorption of Ref. 



3l| . Next we introduce the dispersion and interaction amplitude 



for the BEC-impurity polaron and obtain results for the response to Bragg spectroscopy. 
The result is expressed as a double integral over highly oscillatory integrands that has to be 
calculated numerically and which can give numerical problems. To overcome these problems 
we will then rewrite our expressions in another representation which is suitable for numerical 
work. We then write down the f-sum rule to check our results and find that at low temper- 
ature there is a relatively large recoil contribution for w — i- which numerically has to be 
added separately in the sum rule. This behavior was also found in the optical absorption 
of the solid state Frohlich polaron and in that case there is a relation with the effective 
mass of the polaron j4o|. Recently this sum rule was used to experimentally determine the 
effective mass of the solid state Frohlich polaron, see Refs. 41| and 421]. Next, following De 



Filippis et al 36|, we introduce a lifetime r and examine the dependence of the theory on 
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this lifetime. The results are analyzed for the experimental setup consisting of a Lithium 
impurity in a sodium condensate which is the same system that was examined in Ref. 20 1. 



II. APPLICATION OF THE MORI FORMALISM TO THE DENSITY-DENSITY 
CORRELATION 



The Mori-Zwanzig projection operator technique (see 43|] and 4J]) is a well established 
technique to calculate correlation functions and has been applied before to the Frohlich 
polaron Hamiltonian to calculate the optical absorption by Peeters and Devreese 45| (for a 



detailed description we refer to [46]) . This technique can also be applied to the calculation of 
the density- density correlation function. There is however one obstacle that we encountered 
when defining the Mori projection operator. The reason is that the common definition of 
this operator is ill-defined since it leads to a vanishing density- density commutator. This 
problem was already noticed by Ichiyanagi j47| and he suggested the following definition for 
the action of the Mori projection operator P on an arbitrary operator A: 



PA 



(8) 



In our calculations we have followed this suggestion and found that the density response 
function can be written as: 



(9) 



with: 



O 



Pk^ pI 



Lpk^Pk 
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S f cj, A; 



^ Jo 



dt (1 - e^-*) ( (t) , 4 



(10) 
(11) 
(12) 



Where L is the Liouville operator and S yu, kj is known as the memory function or self 
energy. Furthermore we also introduced the operator Bj: (t) defined as: 



% = QLp^ 



fc' 



(13) 
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with Q the complementary Mori projection operator: Q = 1 — P . The time dependence of 
{t) is governed by a new Liouville operator defined by C =QLQ. It is important to note 
here that until this point no approximations were made and the result for the linear 
density response function is exact. Also, to derive the formula we have not introduced a 
specific system and so it is a general result for the density response function and can also be 
applied to calculate the response of other systems to probes with an arbitrary energy and 
momentum exchange. For example the final state of neutron scattering is also determined 
by the density reponse function 48|. We now consider the Frohlich polaron Hamiltonian 



and follow an approach similar to that used in Ref . |45| . This consists in replacing the new 
Liouville operator C by the sum of the Liouville operator of free Bogoliubov excitations Lboq 
and the Liouville operator of the Feynman modelsystem Lp: C -^Lsog + Lp. Within this 
approximation we can calculate the different quantities in the expression ([9]) which gives 
= and: 



I — iV, 
mj 



(14) 



where is the number of impurities. The memory function is then given by: 



S ( w, A; 
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Im < [e'"''* + 2 cos (ugt) n^] exp 



k + q) D{t) 



(15) 



with the Bose-Einstein distrubution, i.e. the number of Bogoliubov excitations with 
frequency uj^ and D {t) describes the Feynman model system and is given by: 



+ 



2/3 (m + M) 2 (m + M) 
hM 



t 



2mn{m + M) 



1 - exp (int) + 4 sin^ [ — \n{^l) 



(16) 



with M and Q the variational parameters. The values of these parameters are deduced from 
the minimalisation of the free energy 20 1. 



III. LINK WITH THE OPTICAL ABSORPTION OF THE SOLID STATE 
FROHLICH POLARON 

A well-known result from linear response theory and the Kubo formalism is that the 
optical absorption of a system can be expressed as a current-current correlation function 
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49| . Since the current and the density are related to each other through the continuity 
equation it is possible to calculate the current-current correlation function from the density- 
density correlation function. For the process of optical absorption the photon momentum 
is negligible in comparison to the other momenta involved in the polaron problem, which is 
not the case for Bragg spectroscopy. To obtain the optical absorption we have to take the 
limit A; — )■ which gives for the memory function: 

X Im { [e^"^'* + 2 cos {u^t) exp [-q^D {t)] ] . (17) 



This is just the memory function S {u) calculated in Ref. |45| multiplied by uj: 



S (cj, 0) = (w) . 



(18) 



For the solid state Frohlich polaron it is the current-current correlation that is needed for the 
optical absorption. By applying two partial integrations and making use of the continuity 
equation one can find a simple relation between the two: 



die 



iujt 



k.Jt {t) , k.j 



(19) 



We now make use of the standard Kubo formula for the optical conductivity (see e.g. 49|): 
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Re [a (w)] = Re 
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mi a; — E (cj 
this is the result that was found in Ref. 



(20) 



rederived in Ref. 



31| from the Feynman path integral formalism and 



45| using Mori's formalism. 



IV. RESULTS FOR THE BEC-IMPURITY POLARON 

In this section we will introduce the dispersion and the interaction amplitude for the 
BEC-impurity system and obtain an expression for the response of the system to a Bragg 
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pulse. From expression (|5]) it follows that the response of the system is characterised by the 
imaginary part of the density response function which can be written as: 



Im 



A-2 
-—N- 
mi 



Im 



S ( A; 



S ( A; 



+ Im 



2 • 



(21) 



This particular form for the response function is general for the Frohlich polaron and was 
first introduced in 31j. The memory function for the BEC- impurity system is given by (with 



polaronic units i.e. h = mi = ^ = 1): 



f^y r,,^ r.t(i-e-) 



X Im 



A^47r \ mB 

^iuj^t _^ 2 cos (u^t) n^] 



[2kqD {t)f 

exp [-{k + q^Dit)] {[1 + 2kqD {t)f + 1})] . 



(exp [- {k - qf D {t)] { [1 - 2kqD {t)f + l} 



(22) 



The double integral in expression fl22l) could not be calculated analytically and it was com- 
puted numerically. 



V. SECOND REPRESENTATION FOR THE IMAGINARY PART OF THE MEM- 
ORY FUNCTION 



Due to the highly oscillatory behavior of the integrand it turns out that expression ( 122|) 
is not of a form efficient for numerical calculation of the imaginary part of the memory 
function. Therefore we develop an alternative representation for the memory function which 
allows to perform the time integration analytically. This provides a representation for the 
imaginary part that is well-suited for numerical calculations. The real part of the memory 
function turns out to be far more involved but it seems that in this case we can use expression 
(1221) . We will follow a similar approach as was used in Ref. jsil and also in Ref. 45 1. We 
start by rewriting ( fT6i) as: 



D{t) 



h 



2/3 (m + M) 2 (m + M) 



t 



+ 



hM 



2mn{m + M) 



coth 



[l + n (f2)) exp {iVtt) - n (f2) exp {-iVtt) 



(23) 
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After two Taylor expansions for the two exponentials and writing the term with as a 
Gaussian integral we obtain: 



exp \-[k + q) D{t) 



/3 



27r (m + M) 



1/2 



1 1 
n! n'\ 



n,n' 



hM 



2mn{m + M) 



[l + nm 



hM 



2mfi(m + M) 



n{VL) 



2(n+n') /""^ 
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X exp < it 



h[k + q] p 



k + q 



2 (m + M) m + M 



+ {n- n') Vt 



(24) 



Furthermore the memory function can be written as: 



dte'^' / d\{Bj;{t-ih\)Bi 
X Jo Jo ^ 

k.q) /-oo /./? 



Nm 



qj^O 



- (fc + g) D{-t + ih\) 



X exp 

Now the time integration is straightforward 
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k.q 



2(n+n') /"^ 

X U; + g / dpe 2(M+m) 



1-e- 



k^ 



with the following notation: 

a' m 

B{(3,n,n') 



hM 



2mn{m + M) 
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h^VL 



-[a'il + nmYla'nin)]- 



h\k + q 



± 



P 



k + q 



+ {n- n') VL. 



1 _ l + n^ 

hvj^ oj — Vj^ + ie 
(26) 



(27) 
(28) 

(29) 



^ ' 2 (m + M) m + M 

The memory function is now split in an imaginary and a real part with the formula of 
Plemelj. The integral over p in the imaginary part of (|26|) can be performed easily by using 
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the delta function. Taking into account the dispersion and the interaction amphtude, this 
yields (with polaronic units): 



Im 



1 a ( niB + 1 



oo 

X / dqq^ 
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-2qkx)a^(l3) 



1 



g2 + 2 



{qx)' 



X 



^ B {P, n, n) {k^ + + 2qkx) 



(n+n')-l/2 



(1 - e-^'^) 



rig-exp 



/3 (M + 1) A- ju) 
2 (A;2 + g2 + 2gfcx) 



+ (1 + n^) exp 



/3 (M + 1) (u;) 
2 (A;2 + g2 + 2gA;x) 



(30) 



with: 



As shown in Ref. 



(k + q 



+ (n — n') Q — uj. 



45| this expression has a very natural interpretation: every term in 
the double sum corresponds to a well-defined physical process. The (?T,,n')th term in the 
summation represents a scattering process of the polaron where in the initial state the 
polaron is in the n'th internal Franck-Condon state of the Feynman polaron model system 
and is then scattered to the nth internal Franck-Condon state with the absorption of an 
energy huj and a momentum k through the Bragg scattering and with the emission (or 
absorption) of a Bogoliubov excitation with energy hWfj. 



VI. SUM RULE 



As a consistency check of our results we can apply the f-sum rule to our results, this can 
be written as: 



/■oo 

/ dujuj Im 




/o 





2 mi 



(31) 



At small UJ and as temperature tends to zero a delta peak appears in the integrand for 
which, in numerical calculations, we have to add the contribution separately. This feature 
was found for the solid state Frohlich polaron, where the contribution of the delta peak can 



be used to determine the effective mass of the polaron 



40| . The delta peak can be revealed 



if we look at the w — )■ limits of the imaginary and real part of expression (!22|) for the 
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memory function: 



lim Im 
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lim 



S yu, k 



0; 

R [a,k] , 



(32) 
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with: 



47r \ niB 
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X Im 



[e*'^9* + 2 cos (u^t) Uq 



[2kqD {t)f 

exp [-{k + qYD{t)] {[1 + 2kqD {t)f + 1})] . 



(exp [- {k - g)' D (t)] {[1 - 2kqD (t)]' + l} 



If we now look at the w — )■ limit of the density response function we find: 

lim Im \x (uj, k)] = ^ 7r5 (u^) , 

w^-o L V /J mj{l — R{ajj 

where we used the following representation for the delta function: 

1 1 



Six) 
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n^o nvr 1 + (x/n) 
So at temperature zero the sum rule becomes: 
TT k'^N 



mi [l — R[a,k 
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duu Im 
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(34) 



(35) 



(36) 



(37) 



where e is a positive infinitesimal. Care has to be taken at finite temperature since the delta 
peak will then broaden and start to overlap with other contributions of the spectrum. At 
low enough temperatures we will see that our results agree well with fl37|) . 



VII. EXTENDED MEMORY FUNCTION FORMALISM 



It was noticed in Ref. 50| that the Feynman polaron model system does not satisfy the 
sum rule for the density- density correlation. This limitation can be removed by introducing 
a finite lifetime r for the eigenstates of the Feynman model system which mimics scattering 
events with the bosonic bath, which can be done by replacing exp {iVtt) by (1 + it/r)^'^ in 
(fT6|) . This technique is called the extended memory function formalism and was applied in 
Ref. jsGi] for the solid state Frohlich polaron. It was shown that at intermediate coupling 
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this gives a better agreement between the theory and the Diagrammatic Quantum Monte- 
Carlo numerical techniques. This extension also solves the problem of the inconsistency of 
the linewidths at strong coupling with the uncertainty relation for the solid state Frohlich 



polaron, which was first mentioned in Ref. 



3l| . It seems that the second representation 



for the imaginary part of the memory function (1301) is not suitable for the introduction of a 
lifetime. To obtain a qualitatively picture of the dependence of the results on the lifetime we 
can work with the real and imaginary part of the memory function as expressed in equation 

m. 



VIII. RESULTS AND DISCUSSION 

The results we obtain in this section are for a system consisting of a sodium condensate 
with a single lithium impurity, this means that we use for the bosonic mass: ms = 3.8221, 
which is in polaronic units all the results in this section. 

We begin by looking at the weak coupling regime. In figure [2] the Bragg response (12 ip 
is presented as a function of the transferred energy hw for different temperatures and for 
a momentum exchange k = 1. At low temperature (/3 = 100) we clearly see a peak that 
represents the weak coupling scattering process and can be understood as the emission of 
Bogoliubov excitations. Also the contribution of the temperature broadened delta peak at 



low uj is seen. This is the anomalous Drude peak (see Ref. [5l|)- If we look at higher 
temperatures the zero temperature delta peak broadens and there is a larger overlap with 
the weak coupling scattering peak. At /3 = 10 this peak has become indistinguishable from 
the anomalous Drude peak. For /3 = 100 a distinction between the two contributions can 
still be made and the f-sum rule can be applied, as will be done below. 

The dependence of the spectrum on the exchanged momentum k is presented in figure [3] 
for a = 0.1 and /3 = 100. For larger momentum exchange the scattering peak is shifted to 
higher frequencies and a broadening is observed. The inset of figure |3] shows the frequencies 
i^max at which the maximum of the peak occurs as a function of the exchanged momentum 
k together with a least square fit to the Bogoliubov spectrum: 

u = (38) 
2m 

where m is determined as a fitting parameter to be: m = 3.9534; which is in good agree- 
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FIG. 2: The Bragg response of the polaron system consisting of an impurity in a BEG polaron in 
the weak coupling regime (a = 0.1) with momentum exchange k = 1 for different temperatures. 




FIG. 3: The Bragg response of the polaron system consisting of an impurity in a BEG in the weak 
couphng regime (q = 0.1) at temperature (3 = 100 for different momentum exchanges. The inset 
shows the frequency og the maximum as a function of the exchanged momentum k. The curve 
gives a least square fit of the Bogoliubov spectrum (j38p . with fitting parameter m = 3.9534. 



ment with the bosonic mass of the condensate (m^ = 3.8221). This shift according to the 
Bogoliubov dispersion is plausible since the peak corresponds to the emission of Bogoliubov 
excitations. 

In figure H] the high frequency tail of the Bragg spectrum is shown for different coupling 
strengths a at /3 = 100 and k = 1. In the strong coupling regime a resonance is seen 
which is absent in the weak coupling regime. This feature is well-known from the solid state 
Frohlich polaron and corresponds to a transition to the Relaxed Excited State (RES); it 
was first proposed in [52|. This resonance appears at a frequency ujies such that ouj^^s ~ 



Re 



S [ujREs^k 



<^ 1. It is clear 



with the supplementary condition Im 
from f l2ip that these conditions cause a peak in the spectrum. This resonance corresponds 
to a transition from the polaron ground state to an excited state in the polaronic self- 
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FIG. 4: The Bragg response of the polaron system consisting of an impurity in a BEC at different 
values for the couphng parameter a at temperature /3 = 100 and a momentum exchange k = 1. In 
the strong coupling regime the imprint of the Relaxed Excited State resonance appears. The inset 
shows the first three Franck-Condon peaks at a = 8. 



trapping potential which has been relaxed consistent with the new excited wave function of 
the impurity. The coupling strength where the relaxed excited state appears in the Bragg 



spectrum is slightly below a = 4. This is in agreement with the prediction in [20|] that 
for a BEC-impurity the transition between the weak and the strong coupling regime occurs 
around a ^ 3 for /3 — )■ oo. In the strong coupling regime other peaks are present which 
are indicated for a = 8 in the inset of figure HI These are the Franck-Condon (FC) peaks 
and represent a transition to the RES together with the emission of Bogoliubov excitations. 
They only appear in the strong coupling regime which was also observed in the case of the 
acoustic polaron |53l |. 

The dependence of the RES peak on the exchanged momentum is depicted in figure 
[5l The inset shows the frequency of the maximum of the RES peaks as a function of the 
exchanged momentum q together with a least square fit to a quadratic dispersion: 

- = - + ^. (39) 

where i' and m* are the fitting parameters. This suggests that the RES is characterized by 
a transition frequency u and an effective mass m*. 

As a consistency test it was checked whether the above results satisfy the f-sum rule (|37j) . 
Filling out the expression for the imaginary part of the density response function ( 12T|) and 
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FIG. 5: The Realaxed Excited State peak in the Bragg response of the polaron system consisting 
of an impurity in a BEG in the strong coupling regime (a = 8) at temperature (3 = 100 for 
different momentum exchanges. The inset shows the frequency of the maximum as a function of 
the exchanged momentum k at (3 = 100. The curve is a least square fit to the quadratic dispersion 
(j39p . with fitting parameters u = 40.76 and m* = 3.84. 



k 



a 



1 



1 1.5442 1.6845 1.5532 
3 1.5482 1.6419 1.5355 
5 1.5506 1.5555 1.4979 

TABLE I: Results of a numerical calculation of the left hand side of expression (j40p with a cut-off 
coc = 500 introduced for the cj-integral. 

dividing by the common factor the f-sum rule takes the form: 



TT 



2{1- R[a,k 



Im 



dujuj- 



Re 



+ Im 



TT 

2' 



(40) 



S yco, k 

It is impossible to integrate numerically to infinity and for this reason a cut-off, Uc = 500, is 
used. A calculation of the left hand side of f HOj) results in the values in table [T] for different 
values of a and k at (3 = 100. These results should be compared with the rigourous value 
7t/2 = 1.5708... which gives a good agreement with small deviations. 

We now give a qualitative analysis of the dependence of the results on the lifetime r, 
which was introduced in section IVIII within the extended memory function formalism. In 
figure E] the Relaxed Excited State peak is presented at different values for r in the strong 
coupling regime (a = 8). It is observed that the inclusion of a lifetime parameter of the 
order of the polaronic time unit results in a broadening of the peak, where a smaller lifetime 
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FIG. 6: The Relaxed Excited State peak in the Bragg response of the polaron system consisting of 
an impurity in a BEG in the strong couphng regime (a = 8) at temperature (3 = 100 with different 
hfetimes r. 

corresponds to a broader peak. When a lifetime of the order of a hundredth of the polaronic 
time unit is introduced the peak can not be distinguished any more. 

IX. CONCLUSIONS 

We have derived a general formula for the density response function as a function of the 
transferred energy and momentum with the Mori-Zwanzig projection operator formalism. 
This provides a general result that can be used for Bragg spectroscopy but also for other 
probes that exhibit an arbitrary energy and momentum exchange as for example neutron 



scattering where the output is also determined by the density response function 



48|. This 



is applied to the Frohlich polaron Hamiltonian for which the well-known results from Ref. 



3l| for the optical absorption are found. We then extend the analysis to Bragg scattering 
of impurity polarons in a Bose condensed gas, where the Bogoliubov excitations play the 
role of the phonons in the polaron formation. The f-sum rule is checked. 

To analyze the results we introduced the specific system of a lithium impurity in a sodium 
condensate and calculated the spectra in the different coupling regimes and for different 
momentum exchanges and temperatures. It is seen that these spectra possess similar features 
as also found in the optical absorption of the solid state Frohlich polaron. Furthermore it 
was shown that the weak coupling scattering peaks follow the Bogoliubov spectrum as a 
function of the exchanged momentum. In the strong coupling regime the Relaxed Excited 
State emerges, and we derive the transition frequency and the effective mass associated with 
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the Relaxed Excited State. This is of importance for the comparison with Diagrammatic 
Quantum Monte-Carlo numerical techniques, which in the case of the optical absorption of 
the solid state Prohlich polaron has led to new results concerning the linewidth and oscillator 
strength of the Relaxed Excited State and Franck-Condon transitions. The Franck-Condon 
peaks were also observed. 

Our results were tested using the f-sum rule which resulted in a good agreement with 
small deviations. 

The influence of the introduction of a lifetime within the extended memory function for- 
malism was also qualitatively investigated and it was shown that this results in a broadening 
of the RES peak. 
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